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I Abstract. We present the results of a BeppoSAX obser- 
i vation of the nearby spiral galaxy NGC 4945 in the 0.1- 
I 200 keV energy band. The nuclear X-ray emission emerges 
. above ~7 keV, through an absorber with column density 
] Nh ~ a few 10^'' cm~^. Its remarkable variability (dou- 

■ bling/halving timescale ^lO** s) is not associated with 
I any appreciable spectral changes, ruling out variations of 
I the interposed absorber or changes in the primary contin- 
' uum shape. The intense iron Kq, fluorescent emission line 
I. is likely to be produced by the same nuclear absorbing 
' matte. Our estimate of the intrinsic nuclear power sug- 
. gests that the active nucleus dominates the bolometric en- 
I ergy output ofNGC 4945, if its nucleus has a typical quasar 
. Lx/Lboi ratio. This supports the idea that X-rays are the 

■ best energy band to search for absorbed AGN. The 0.1- 
5 keV emission is extended along the plane of the galaxy, 
and most likely due to a population of unresolved bina- 

. ries, as believed to happen in several early-type galaxies. 
I Alternatively, hot gas associated with a starburst-driven 
. superwind outflow could substantially contribute to the 
extended emission above 1 keV. 



Key words: Galaxies: active - Galaxies: individual: 
NGC 4945 - Galaxies: nuclei - Galaxies: Seyfert - X-rays: 
galaxies 

1. Introduction 

NGC 4945 is a nearby (z=0.0019 or 3.7 Mpc; Mauers- 
berger et al. 1996), almost edge-on (inclination angle 
~ 80°) spiral galaxy, which exhibits a prominent dust lane 
crossing its plane. It is the third brightest extragalactic 
source in the IRAS point source catalog, and most of its 
infrared emission is concentrated in a compact nuclear re- 
gion (Rice et al. 1988; Brock et al. 1988). It shows both 
starburst emission (Heckman et al. 1990; Koornneff 1993) 
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and an H2O megamaser (Dos Santos & Lepine 1979). 
A region of extended optical-line emitting gas protrudes 
from the nucleus along the galaxy minor axis (Nakai 1989; 
Heckman et al. 1990). It is associated with a cavity open 
by a starburst superwind (Moorwood et al. 1996). X-ray 
observations with Ginga (Iwasawa et al. 1993) and ASCA 
(Iwasawa 1997) unveiled a highly obscured, strongly vari- 
able X-ray source. Actually NGC 4945 turned out to be 
the brightest Seyfert 2 and the second brightest radio- 
quiet Active Galactic Nucleus (AGN) after NGC 4151 of 
the 100 keV sky (Done et al. 1996). 

NGC 4945 is one of the best examples of an active 
galaxy with a composite nature (starburst plus AGN) . At 
wavelengths shorter than 1 keV, all its observational prop- 
erties can be accounted for by starburst activity alone, 
despite its classification as a Seyfert in the Veron-Cetty 
and Veron (1989) catalog. For instance, a mass-to-light 
ratio of 0.18 is consistent with the parameter space solely 
occupied by starbursters (Oliva et al. 1999). Also in the 
Genzel et al. (1998) ISO diagnostic planes NGC 4945 is 
located in the regions occupied by starbursts. However, 
the presence of a strong continuum radio source and stud- 
ies of the off-nuclear optical spectra have suggested the 
presence of a LINER (Whiteoak & Gardner 1979; Moor- 
wood et al. 1996). The X-ray observations provided the 
final proof for the presence of an active nucleus. This con- 
firms that hard X-ray emission can be the most efficient 
wavelength to identify the nature of obscured AGN (see 
e.g. the discussion in Vignati et al. 1999 for the case of 
NGC 6240). 

Previous observations in the soft (0.1-2 keV) and in- 
termediate (2-10 keV) X-rays showed a complexity, which 
is far from being completely resolved. The 3-10 keV Ginga 
spectrum (Iwasawa et al. 1993) could be fitted with a 
simple power-law with a photon index of about 1.7, in- 
terpreted as scattering of the primary nuclear contin- 
uum. The same description fits also the non-simultaneous 
Ginga, ASCA and OSSE spectra (Done et al. 1996). A 
huge (Equivalent Width, EW, ^^1-1. 5 keV) iron line has 
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been observed as well, whose eentroid energy is eonsistent 
with Kq, fluorescence from neutral or mildly ionized iron. 
The ASCA image is clearly extended, but the irregular, 
broad shape of the instrumental Point Spread Function 
(PSF) has prevented any detailed studies. The ROSAT 
PSPC and HRI images of NGC 4945 revealed a complex 
pattern. Brandt et al. (1996) detected at least five discrete 
sources, one of them variable by about one order of magni- 
tude on timescales of hundreds days. A diffuse emission up 
to a distance of about 5' from the nucleus, and elongated 
along the plane of the galaxy, was measured as well. 

The scientific payload on board BeppoSAX (Boella et 
al. 1997a) is particularly well suited to study this source 
and resolve at least part of this complexity. It combines 
the widest X-ray broadband coverage, the highest sensi- 
tivity above 10 keV and the sharpest 2-10 keV PSF ever 
flown before the advent of Chandra. We report in this pa- 
per the results of a BeppoSAX observation performed on 
July 1999. At the distance of NGC 4945, 1' corresponds 
to about 1.1 kpc. 

In this paper: energies are quoted in the source 
rest frame; uncertainties on the spectral parameters 
are quoted at the 90% level for one interesting pa- 
rameter (Ax^ = 2.71); the cosmology used assumes 
Ho = 50 km s~^ Mpc~^ and qo = 0.5; J2000 coordinates 
are used, unless otherwise specified. The Galactic col- 
umn density along the line of sight is assumed to be 
1.57 X 10^1 cm-2 (Hciles & Cleary 1979). The errors on 
the energies measured by the MEGS take into account a 
systematics of 0.8% at 6 keV (Guainazzi & Molendi, 1999). 
XsPEC 10 was used for spectral analysis. 

The paper is organized as follows. In Sect. 2 the data 
reduction procedures are described. In Sect. 3 we will deal 
with the analysis of the LEGS and MEGS images. We 
will report the detection in the intermediate X-rays of two 
discrete sources in the SE outskirts of the NGG 4945 spa- 
tial profile (Sect. 3.1), and study the residual unresolved 
diffuse emission (Sect. 3.3). We present also a reanaly- 
sis of ROSAT/PSPG archival data of the NGG 4945 field 
(Sect. 3.2). The properties of the hard {i.e.: >10 keV) X- 
ray emission will be studied in Sect. 4.1. In Sect. 4.2 we 
will describe the broadband X-ray spectrum with standard 
fitting techniques. We will discuss our findings in Sect. 5 
and summarize them in Sect. 6. 

2. Observation and Data Reduction 

NGG 4945 was observed by BeppoSAX from July 1 04:43 

UTG to July 3 1999 09:18 UTC. The scientific payload 
onboard BeppoSAX comprises two gas scintillation pro- 
portional counters with imaging capabilities: the Low En- 
ergy Goncentrator Spectrometer (LEGS; Parmar et al. 
1997), and the Medium Energy Goncentrator Spectrome- 
ter (MEGS; Boella et al. 1997b). The LEGS has a nomi- 
nal bandpass of 0.1-10 keV, with an energy resolution of 
-4% at 1 keV and -8% at 6 keV. The field of view has 



a diameter of ~37'. The angular resolution at 2 keV is 
2'.1 FuU Width Half Maximum (FWHM), but degrades 
to 9'. 7 below the Garbon edge {i.e., <, 0.4 keV). The 
MEGS has similar energy resolution, but a wider field of 
view (57' diameter), a narrower sensitive bandpass (1.8- 
10.5 keV) and about a factor of two higher effective area 
in the overlapping energy interval, after the failure of one 
of the three original units in May 1997. The 80% power 
radius of the PSF is comprised between 2'. 5 and 2'. 75. 
The BeppoSAX payload also includes two collimated in- 
struments, mounted on a rocking system to achieve a con- 
tinuous monitoring of the background, with a duty-cycle 
of 96 s. The High Pressure Gas Scintillator Proportional 
Gounter was switched off during the NGG 4945 observa- 
tion. The Phoswitch Detector System (PDS; Prontera et 
al. 1997) possesses an unprecedented sensitivity in its 13- 
200 keV nominal bandpass. 

The data reduction followed standard procedures, as 
described e.g. in Guainazzi ct al. (1999b). In particular, 
good time intervals for the extraction of scientific products 
were selected only when the star tracker aligned with the 
pointing direction was operative. This condition allows the 
best attitude reconstruction. The absolute pointing accu- 
racy under these conditions is better than 1', the relative 
better than 30" for sources as bright as NGG 4945. The 
resulting net exposure times are about 38.8 ks, 86.7 ks 
and 87.6 ks in the LEGS, MEGS and PDS, respectively. 
PDS data have been reduced using fixed Rise Time thresh- 
olds to reject the particle background, as appropriate for 
sources brighter then 0.5 s~^. The PDS points continu- 
ously at the nominal target with only two of its four units, 
the others monitoring the background. 

MEGS background spectra were extracted from blank 
deep field exposures, accumulated by the BeppoSAX Sci- 
ence Data Genter (SDG) in the first three years of the 
mission. The same method applied to the LEGS spec- 
tra yields systematically negative coimts below the car- 
bon edge. We have therefore extracted LEGS background 
spectra from two semi-annuli in the LEGS field of view, 
and siutably renormalized them to the source extraction 
region, as described in Parmar et al. (1999). The back- 
ground subtraction in PDS light curves and spectra has 
been performed by plain subtraction of the "off-source" 
from the "on-source" products. The systematic uncertain- 
ties of this method are lower than ~ 0.03 s~^ in the 
full PDS sensitive energy bandpass (Guainazzi & Mat- 
teuzzi 1997). Spectra and light curves of NGC 4945, un- 
less otherwise specified, have been extracted from circular 
regions of 8' radius around the best-fit eentroid of the 
2-10 keV image. The position of the eentroid is slightly 
dependent on energy. However, the dynamical range of 
these fluctuations is of the order of 20", well within the 
accuracy of the positional reconstruction. The spectra 
were rebinned in order to: a) oversample the FWHM 
of the energy resolution by a factor at least 3; b) have 
at least 20 counts per spectral channel, to ensure the 
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Fig. 1. Iso- intensity contours for the MECS 2-10 keV im- 
age (solid lines) superimposed to the ESO DSS optical 
image. The contours correspond to steps of la of the fo- 
cal background fluctuations starting from Aa. A Gaussian 
smoothing with width equal to 2.5 pixels has been applied 
to the MECS image for display purposes only. The best-fit 
centroid positions of the serendipitous sources SI and 82 
are marked 

applicability of statistics. The background-subtracted 
count rates are: (2.51 ± 0.10) x 10^^ s'^ (0.1-4 keV); 
(6.30 ± 0.10) X 10-2 s-i (1.8-10.5 keV); 2.72 ± 0.04 s^i 
(13-200 keV) in the LEGS, MEGS and PDS, respectively 
(in brackets the energy intervals where each instrument 
is currently well calibrated). The September 1997 calibra- 
tion release is adopted throughout this paper. No known 
bright sources are present in the PDS field of view apart 
from NGG 4945 itself, and the probability of a serendipi- 
tous source with a flux equal or larger than NGG 4945 is 
^ 2.5 X 10-4, if the 2-10 keV ASGA LogN-LogS (Gagnoni 
et al. 1998) is adopted and the lowest extrapolation of the 
unabsorbed 2-10 keV flux is assumed (see Sect. 4.1). 

3. Image Analysis in the 0.1-10 keV Energy Band 

3.1. The serendipitous sources 

The rather large visual extent of NGG 4945 (3'x20') 
is well encompassed within the fields of view of both 
the LEGS and the MEGS. It is therefore possible with 
BeppoSAX to compare directly the shapes of the opti- 
cal and of the X-ray images. In Fig. ^ we superpose the 
2-10 keV (MEGS) iso-intensity contours with the ESO 
Digitalized Sky Survey image. The MEGS contours are 
well aligned with the main plane of the host galaxy. Ac- 
tually, part of the X-ray extension is likely to be due 
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Fig. 2. Serendipitous source SI spectrum [upper panel) 
and residuals in units of standard deviations (lower panel) 
if a power-law model absorbed by a column density 
Nh = Nn.Gai is applied 



to the presence of discrete sources. Apart from the nu- 
cleus, one source is detected at a signal-to-noise ratio 
>5 (SI hereinafter). Its coordinates are a = 13*^05^6. 6"; 
5 = -49°33'34", and it is therefore located 6.2' SW from 
the nucleus. Assuming Nh = Nn^cai (the 90% upper limit 
on Nh if left free in the fit is 1.1 x lO^"^ cm-^), its spec- 
trum (extracted from a 2' circular region) can be mod- 
eled with a power-law with T = 1.7±j]:| (x^ = 4.4/9 dof; 
see Fig. ^) or a thermal bremsstrahlung with kT > 4 keV 
(x^ — 3.8/9 dof). The net background subtracted count 
rate is (1.7 ± 0.2) x 10-^ s-^ in the 2-10 keV band, cor- 
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responding to a flux of ~ 1.7 x 10 
rest-frame luminosity of ~ 3 x 10"^^ erg s 
of NGG 4945. The search for a periodic modulation has 
produced no significant peak, hardly surprising, since only 
147 photons are counted from SI. Another detection at the 
4ct level [count rate (7.3 ± 1.7) x IQ-^ s-^] corresponds to 
the following position: a = 13''04'"41.9^ S = -49°34'22" 
(9' SW from the center). The faintness of this source, and 
its location below the MEGS strongback rib prevent us 
from performing any reliable flux determination or spec- 
tral analysis. 

In Fig |[ we show the NGG 4945 field imaged by the 
ROSAT/PSPG on August 1992 and July 1993, for a to- 
tal net exposure time of 14.2 and 9.0 ks, respectively. The 
BeppoSAX sources coincide with two of those detected by 
ROSAT within the 1' positional accuracy. The extrapola- 
tion of the SI MEGS best-fit into the PSPG band predicts 
48 ±6 photons, reasonably close to the observed 33 PSPG 
background-subtracted source counts. Despite the rela- 
tively high Galactic latitude of NGG 4945 (b = 13.340), 
the possibility that either or both SI and S2 is a Galac- 
tic foreground object cannot be ruled out by our data 
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Fig. 3. Adaptively smoothed image of the ROSAT/PSPC 
NGC 4945 field. The circles indicate the position of SI 
{left) and S2 {right) and have a radius of 1', corresponding 
to the typical MEGS error circle. The white crosses mark 
the position of the sources hsted in Tab. 



alone. No known source is included in the Simbad or Ned 
databases within 2' from the best-fit MEGS positions. If 
Nh = 10^° cm~^ (appropriate for a Galactic source at the 
distance of approximately 0.5 kpc), the extrapolation of 
the MEGS best-fit would exceed the ROSAT flux by a 
factor of four, which could be accounted for by intrinsic 
variability. Analogously the extrapolation of the MEGS 
spectrum remains consistent with the ROSAT detection 
within a factor of a few for absorbing column density 
<! 10^^ cm"^. It is therefore also possible that SI is a 
background object, seen through a moderate column den- 
sity of the galaxy. 



3.2. The ROSAT nuclear sources 

The ROSAT PSPG detects four sources within the inner- 
most NGG 4945 5' (see Fig. |). Their coordinates and 
count rates are listed in Tab. l|. All of them exhibit a 
certain amount of variability. It is interesting to notice 
that none of them strictly coincides with the optical nu- 
cleus of the galaxy, the closest being located 0'.6 SE. 
An extended diffuse emission along the galactic plane 
is also present, whose integrated count rate is about 3- 
4x10"^ {i.e., about the same as the sum of the de- 
tected discrete sources). No reliable spectral information 
can be obtained from any of the above sources. If their 
spectrum is a power-law with F = 2 (1), they may con- 
tribute -15% (50%) of the observed 2-10 keV MEGS flux. 
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"0.1-2.4 keV PSPC background subtracted count rate 

''August 1992 observation 

'^July 1993 observation 

''distance from the galaxy optical nucleus 

■^not detected 

Table 1. Properties of the sources detected in the 
ROSAT/PSPG NGG 4945 field within 5' from the opti- 
cal nucleus 



3.3. The diffuse emission in BeppoSAX 

In order to provide a quantitative estimate of the in- 
trinsic extension and ellipticity of the BeppoSAX X- 
ray images, we extracted the counts along two perpen- 
dicular strips crossing the optical center of the galaxy 
in the 0.1-2 keV (LEGS) and 2-10 keV (MEGS) en- 
ergy bands, respectively. The strip directions are defined 
by the following pair of coordinates: {a = 13''05™44.3''; 
S = -49°25'22") and (a ^ 13''05"n2.3"; S = -49°30'34") 
for the strip parallel to the plane of the galaxy {i.e.: 
NE to SW); (a = 13'^05"^36.P; (5 = -49°30'26") and 
{a = 13'^05™13.1''; S = -49°26'10") for the strip perpen- 
dicular to the plane of the galaxy {i.e.: NW to SE). The 
width of the strips is 4', comparable with the intrinsic 
width of the instrumental PSF (Parmar et al. 1997; Boella 
et al. 1997b). This choice is actually likely to discard the 
bulk of the photons emitted below the carbon edge, where 
the LEGS detector PSF is the broadest. However, this ef- 
fect is negligible in our case, because the photoelectric 
absorption due to the matter in our Galaxy along the 
line of sight of NGG 4945 is large enough to substan- 
tially suppress the photons in this band. The results are 
shown in Fig. ^ as count rate space density as a func- 
tion of the offset angle towards the galaxy center (note 
that the offset increases for decreasing RA). The MEGS 
image shows clearly a larger extent in the direction par- 
allel to the galaxy plane for all offset angles between -4' 
and 8'. Gomparison at larger angles is made impossible by 
the shadow of the strongback support ring in the MEGS 
field of view. The difference between the profiles along 
the parallel and perpendicular directions is rather sym- 
metric within 4', while an excess tail is evident at larger 
offset radii only for positive offsets (cf. the inset in the 
right panel of Fig. ^) , due to the presence of the discrete 
sources detected by the MEGS. 

In the same Fig. ^ we compare the observed NGG 4945 
profile with point-like sources, whose data have been re- 
duced under the same experimental conditions as the 
NGG 4945 ones. In order to avoid any systematic effects 
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Fig. 4. Count rate profiles along directions parallel {empty circles) and perpendicular {filled squares) to the optical 
plane of NGC 4945 for the LECS 0.1-2 keV {left panel) and the MECS 2-10 keV {right panel) images. The solid lines 
represent the radial profile for the BeppoSAX observations of LMC X-3 and 3C273, respectively. The offset positive 
sense is taken from NE to SW, therefore in the sense of decreasing RA. In the inset of the right panel the quantity A 
is plotted as a function of offset angle, defined as the difference between the observed count rates in the parallel and 
perpendicular directions at the same offset distance. Each tickmark in this inset corresponds to 2' 



due to the dependence of the PSF with energy, we have 
chosen two sources contained in the archive of the Bep- 
poSAX public observations, whose spectrum is similar to 
that observed in NGC 4945. The 0.1-2 keV LECS spec- 
trum of NGC 4945 can be formally well approximated by 
a photoelectrically absorbed power-law, with parameters: 
r = 1.2±°;3; Nh = (1-2±J:|) X 10^1 cm-2. This model is 
consistent, within the statistical uncertainties, with that 
describing the LECS spectrum of the black-hole candi- 
date LMC X-3 in the same energy band (Siddiqui et al. 
1999). The 2-10 keV MECS continuum of NGC 4945 is 
well approximated by a power-law, with F = 1.91±Q'^g. 
The comparison for the MECS is made against the profile 
of the quasar 3C273 (Grandi et al. 1997). The presence of 
a prominent iron line in the NGC 4945 spectrum (Iwasawa 
et al. 1993), which is absent in 3C273, makes a negligible 
difference in this context. Even the profile perpendicular 
to the galaxy plane is broader than the instrumental PSF 
within ±3'-4'. However, caution has to be used in inter- 
preting the last result, as the residual inaccuracies in the 
attitude reconstruction could significantly affect the sig- 
nificance of the last result, whereas no doubt exists about 
the true nature of the extension along the galaxy plane. 

As a further step, we have studied the energy de- 
pendence of the X-ray image profile along the plane of 
the galaxy with the MECS (an analogous study with the 
LECS is hampered by the much poorer statistics and the 
broader instrumental PSF). The results of such a study 
are shown in Fig. |[ They are compared with the width 
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Fig. 5. FWHM of the X-ray image profile along the host 
galaxy plane as a function of energy. The dash-dotted line 
marks the values expected from the on-flight MECS cali- 
bration (Boella et al. 1997b; Chiappetti 1998). In the inset, 
the percentage fractional difference between the observed 
and the expected values as a function of energy is shown. 

expected on the basis of the instrumental PSF, weighted 
on the energy assuming a power-law spectrum of photon 
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Fig. 6. Comparison of the MECS and PDS spectra of 
NGC 1068 {filled squares) and NGC 4945 {empty circles). 
Both spectra have been normaUzed to the Crab spectrum, 
which is a power-law with a photon index ~2. Each data 
point corresponds to a signal-to-noise ratio >10, except 
the PDS data points of NGC 1068, for which the signal- 
to-noise ratio is >3. In the inset a zoom of the 5-8.5 keV 
energy band, with the flux in linear scale 

index T — 1.8 (dash-dotted line in Fig. ||). The X-ray 
MECS image is significantly extended in the whole range 
between 2 and 5 keV, whereas it becomes consistent with 
that expected from a point-like source for E ^ 6 keV. It 
is difficult with our data to assess a systematic trend be- 
tween 2 and 5 keV. If we consider the face-on values, the 
profile width has a peak at ~4 keV. 

4. The X-ray Spectral Energy Distribution 

In Fig. H we compare the >2 keV spectra of NGC 4945 
and NGC 1068, after they have been divided by the spec- 
trum of the Crab Nebula. Since the latter is a power-law 
with photon index ~2, this technique is very similar to a 
vFi, plot. In the intermediate X-ray band, the two spectra 
are similar. A careful inspection suggests that the emis- 
sion line feature around 6 keV is broader and more intense 
in NGC 1068, indicating that more ionization stages con- 
tribute to its profile than in NGC 4945. The most remark- 
able difference is, however, at energies >10 keV. NGC 4945 
is more than one order of magnitude brighter in the PDS, 
suggesting that a further component emerges - actually 
well connected with the higher end of the MECS sensitive 
bandpass. It is straightforward to identify this component 
with the primary nuclear continuum transmitted through 
an absorber with a column density of a few 10^"* cm~^, 
originally discovered by Iwasawa et al. (1993). On the 
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Fig. 7. Light curve in the 13-200 keV energy band (PDS). 
The binning time is 11520 ks, approximately correspond- 
ing to two BeppoSAX orbits 

other hand, in NGC 1068, the archetypical example of 
totally Compton-thick Seyfert 2, the nuclear emission is 
completely suppressed. 

4-.1. The high- energy emission 

As shown in Fig. 0, the PDS full band light curve exhibits 
a remarkable variability. In three events the flux increases 
by about 60% in ~3 x lO'' s, whereas a decrease by a factor 
~30% in <i 10^ s represents the most extreme variability 
episode. We have searched in vain for spectral changes 
associated with these flux variations. In Fig. || we show the 
hardness ratio (HR) between the count rates in the 13-35 
and 35-200 keV energy bands. No clear deviation from 
constancy is observed. If we fit the data points in Fig. |^ 
with a constant, = 18.0/17 degrees of freedom (dof). 
We will therefore hereinafter focus on the time-averaged 
spectrum only. 

An absorbed power-law model yields a marginally ac- 
ceptable fit to the PDS spectrum (x^ = 62.3/47 dof). 
We first modeled the extinction as an exponential func- 
tion of the energy exp"'^"'^^^), where the cross-section 
cr(E) includes both the effects of photoelectric absorp- 
tion (model wabs in XSPEC, with the abundances of Mor- 
rison & McCammon 1983) and of Compton scattering 
(model cabs in XsPEC). The inclusion of a high-energy 
cutoff improves the quality of the fit by Ax^ — —7.7, for 
the decrease of one degree of freedom, which is signifi- 
cant at > 99.2% confidence level according to the F-test. 
The best-fit parameters and results are summarized in 
Tab. 1^. In Fig. |^ we show the iso-x^ confidence contours 
for r versus the cutoff energy Ec. At 90% confidence level 
for two interesting parameters, Ec is comprised between 
70 and 500 keV, the value corresponding to the best-fit 
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Fig. 8. Hardness ratio between the count rates in the 35- 
200 keV and 13-35 keV energy bands, as a function of the 
13-200 keV count rate 



Model 


Nh 
(10^^ cm-2) 


F 


(keV) 


xVdof 


PL" 

CP" 
CP*" 


5.4 ±0.8 
3.9±!]| 
2.2±[];1 


1.4 ±0.3 


110±^g« 

i30±iE;'' 


62.3/47 
54.7/46 
57.7/46 



"Compton absorption after model cabs in Xspec 



''Compton scattering after Matt et al. (1999b; absorption 
model for optically thick Compton scattering). In comparison 
to the same model without high-energy cutoff Ax^ = -6.9 
Table 2. Best-fit parameters for the fits on the PDS spec- 
trum alone. PL = power-law, CP = cutoff power-law. 



being ~110 keV. The addition of a Compton-reflection 
component from a plane-parallel infinite edge-on slab to 
the best-fit model (implying the decrease of a further de- 
gree of freedom to the fit, the relative normalization be- 
tween the reflected and the transmitted component R; 
model pexrav in XsPEC, Magdziarz & Zdziarski 1995) 
does not yield any further improvements to the fit. The 
90% upper limit on R is 0.5. The 20-100 keV observed 
flux is (2.89 ±0.04) x 10"^° erg cm'^ 8"^ The extrap- 
olation of the best flt model (second row of Tab. |^) to 
the 2-10 keV (0.1-200 keV) energy band yields an unab- 
sorbed flux of ~ 3.5 x 10"^ (1.83 x 10"*) erg cm'^ s-\ 
corresponding to a rest-frame luminosity of ~ 5.3 x 10"*^ 
(2.8 X 10^4) erg s-\ 

At the level of 4 x 10^^ cm~^, the matter is opti- 
cally thick to Compton scattering. The cabs model in 
Xspec assumes, on the contrary, the simple Thompson 
cross-section. Neglecting the effects of electron scattering 
in the column may lead to a substantial overestimate of 




o 



Fig. 9. Iso-x^ confidence contours for the spectral index 
versus the high-energy cutoff, for the best-fit model of the 
second row of Tab. 1. Contours are at 67%, 90% {solid 
lines) and 99% (dotted lines) for two interesting parame- 
ters. 

the true column density and, consequently, of the intrin- 
sic nuclear luminosity (Leahy et al. 1989). Matt et al. 
(1999b) have recently developed a self-consistent model 
for X-ray absorption by matter which is optically thick 
to Compton-scattering and has a large covering factor. 
The application of such a model yields about 40% lower 
value for the column density. Moreover, the extrapolated 
intrinsic luminosities differ significantly. In the Matt et 
al. (1999b) model a fraction of the photons are scattered 
into the line of sight, and the nuclear luminosity required 
to explain the observed flux is lower. The extrapolation 
of the best flt model (row 3 in Tab. ^) yields an un- 
absorbed 1-10 keV (0.1-200 keV) flux of ~ 2.5 x IQ-^" 
(1.14 x 10~^) erg cm~^ s~^, corresponding to a rest frame 
luminosity of ~ 3.8 x 10'*^ (1.77 x lO'^^) ej.g g-i^ 
refer to these values hereinafter. 

4-2. The broadband 0.1-200 keV spectrum 

Given the spatial and spectral complexity emerging from 
the above analysis, a global description of the 0.1-200 keV 
spectrum in terms of standard fitting technique at the 
moderate spatial resolution and sensitivity provided by 
the BeppoSAX instruments is hardly more than an aca- 
demic task. The soft and intermediate X-ray emission is 
likely to be produced by the superposition of radiation 
scattered/reprocessed in the nuclear environment (prob- 
ably most contributing to the huge iron line), a popula- 
tion of discrete unresolved sources (binaries, cataclysmic 
variables, supernova remnants) and truly diffuse gas emis- 
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Fig. 10. Light curves in the 0.1-2 keV, 2-4 keV and 4- 

10.5 keV energy bands, and hardness ratios between ad- 
jacent bands. The binning time is 11520 s. HRi is defined 
as the ratio between the counts in the 2-4 keV and 0.1- 
2 keV energy bands. HR2 is defined as the ratio between 
the counts in the 4-10.5 keV and 2-4 keV energy bands. 

sion, each possibly seen through absorbers with different 
density and/or physical conditions. Such complex scenario 
will be fully understood when future observation with arc- 
second spatial resolution and/or much better sensitivity 
will be available. Given all these caveats, in this Section 
we will, however, try to accomplish this task, whose results 
has allowed us to extend previous studies of this kind. 

In Fig. ^ we show the light curves in the 0.1-2 keV, 
2-4 keV and 4-10.5 keV energy bands, along with the 
corresponding HRs between adjacent bands. No significant 
variability of either quantity is observed. This is of course 
not surprising, because the bulk of the emission below 
6 keV is extended on several kpc scale (see Sect. 3.3). Fit 
with a constant line of the HRs yield of 22.3/17 dof and 
17.9/17, for HRi (ee2-4 keV/0.1-2 keV) and HR2 (=10.5- 
4 keV/2-4 keV), respectively. We are therefore justified in 
focusing on the time-averaged spectrum only. 

The broadband spectrum of NGC 4945 was fitted using 
the following general formula: 

F{E) = exp'['^'"'^«'G»'l{exp-['^'"-(^)^«'il A{E) 

+ exp-['^'"'(^)^«-=^l B{E) + G{E) 

The last term in brackets is the nuclear transmitted 
component emerging in the PDS, with normalization N. 
fph(E) and crc(E) are the photoelectric and Compton 
cross-sections, respectively. We will use the cabs Xspec 
implementation for the latter. A(E) and B(E) are two con- 
tinuum models, chosen according to one of either the fol- 
lowing scenarios: 



a) A(E) is a power-law, produced by the scattering of 
the primary continuum by ionized matter, and B(E) 
is the thermal emission from a collisionally excited 
plasma. The index of the former is tied to be the 
same as Ftr. We will use the mekal Xspec imple- 
mentation for the latter throughout this paper. This 
model provides a reasonable description of the spectra 
of some reflection-dominated Seyfert 2 galaxies, like 
Circinus Galaxy, NGC 1068 (Guainazzi et al. 1999a), 
or NGC 6240 (Iwasawa & Comastri 1998; Vignati et al. 
1999). All these objects are characterized by a strong 
circumnuclear starburst and/or are infrared ultralumi- 
nous galaxies 

b) A(E) is a thermal plasma emission and B(E) is either 
a blackbody (bl) or a second thermal plasma (b2), 
which is the template model successfully employed to 
fit the Einstein, ROSAT and ASCA bulge emission of 
early-type galaxies (Bregman et al. 1995; Matsumoto 
et al. 1997; Kim et al. 1996; Irwin & Bregman 1999) 

c) A(E) is a power-law and B(E) is a Compton-reflection, 
both with the photon index tied to be the same as Ftr- 
This scenario follows the idea that the primary nu- 
clear continuum undergoes both scattering by a warm 
plasma and Compton reflection, along two different op- 
tical paths. 

To all models a Gaussian emission line G(E) was 
added. We allowed different absorbing column densities 
for the nuclear (warm-scattered and Compton-reflected 
power-law) and the thermal/non-nuclear components. Fi- 
nally, we included in the models a power-law, with all the 
parameters frozen to the best-fit values obtained fitting 
the SI spectrum alone (this contribution is not explicity 
indicated in the above formula for sake of clarity) . 

The best-fit parameters and results are shown in 
Tab. ||. The fit is actually only marginally acceptable for 
scenarios a) and b). However, no significant improvement 
is obtained if a further component (either a continuum, an 
emission line or an absorption edge) is added to any of the 
above models. No improvement in the quality of the fit in 
scenario a) is obtained if the warm scattered power-law 
spectral index is untied from Ftr- No systematics trend 
can be recognized in the residuals, as shown in Fig. pj] . 
We checked that the rather high is not due to residual 
systematics effects . 

The fit results strongly argue against any Compton 
reflection component to substantially contribute to the 
broad band spectrum. The scenario c), which includes 
it explicitly, yields a significantly worse with a com- 
parable number of degrees of freedom. The addition of 
a Compton reflection component to the other scenarios 
is not statistically justified. In scenario a), its 2-10 keV 
fiux upper limit is 5.9 x 10^^'^ erg cm^^ s~^, against a 
3.1 X 10~^^ erg cm~^ s~^ flux of the warm scattered com- 
ponent in the same energy band. 
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Thermal components Power-law 





Nh,i 


kT 


Z 


kT 


Nh,2 


r 




NH.tr 


Ecutoff 






(10^1 cm-2) 


(keV) 




(keV) 






(%) 


(10^* cm-2) 


(keV) 




a) 


<3.4 


3±§ 






26 ±14 


1.61 ±0.16 


0.4/<l'' 
0.4/<12'' 
.../<12'' 

0.4/<2 


4.5 ±0.7 


190±^r 


158.3/127 


bl)" 


6±3 


8±2 


0.14±8:?t 


o.ii±g:St 




i.59±[]:g 


4.5 ±0.7 


140±ir 


158.3/127 


b2)= 
c) 


7±4 


8±2 


o.i2±g:f^ 


o.26±8:?f 




i.6±8:i 
i.70±8:JJ 


4.4±8:i 
5.i±8:i 


140±i^» 
170±JJ» 


155.3/125 
172.2/128 



"warm scattering fraction 

''upper limit on the Compton reflection relative normalization parameter. It is calculated after adding this component to the 
best-fit model 

°in scenario bl) an optically thin plasma thermal emission (columns 2-4) and a blackbody (column 5) are employed. In scenario 
b2) the blackbody is substituted by an optically thin plasma thermal emission 

Table 3. Best- fit parameters and results for the continuum components in the broadband fits of the BeppoSAX 
NGC 4945 observation 




Energy (keV) Energy (keV) 

Fig. 11. Spectra and residuals in units of standard deviations {upper panels) and corresponding best-fit Spectral 
Energy Distributions {lower panels) for scenarios a) {left) and bl) {right) 
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Model 


E 




I 




(keV) 


(eV) 


(10^^ photons cm"^ s"^) 


a) 


6.46 ± 0.07 


<150 


4.3±?J 


bl) 


6.46 ± 0.07 


<140 


4.2 ± 1.1 


c) 


6.48 ±0.07 


<140 


4.2 ±0.7 



Table 4. Iron line properties for the best-fit models of 
Tab. 3 



The properties of the iron line are not strongly de- 
pendent on the adopted continuum (see Tab. U). It re- 
mains narrow (intrinsic width cr < 150 eV) and consistent 
with Kq fluorescence from neutral or mildly ionized iron 
(Fe<xx). If we extrapolate the best-fit model of raw 3 
Tab. ^, and compare it with the observed line intensity, 
the iron line EW against the transmitted nuclear contin- 
uum alone is ~1.3 keV. 

Assuming the best- fit scenario b), the observed flux 
in the 0.1-2 keV (2-10 keV) energy band is 1.3 x lO'^^ 
(5.4 X 10~^^) erg cm~^ s~^, corresponding to a rest frame 
unabsorbed luminosity of 1.7 x 10"*° (8.4 x lO'*'^) erg s~^. 

5. Discussion 

5.1. The hard X-ray emission 

NGC 4945 is confirmed to be one of the brightest extra- 
galactic objects above 10 keV, where we are seeing the 
primary nuclear continuum, transmitted through an ab- 
sorbing column density of a few 10^^ cm~^. This interpre- 
tation is confirmed by the detection of rapid variability in 
the PDS emission, with a extrapolated doubling/halving 
time scale r ~ 3-5 x lO'' s. Variation of a comparable 
amount were observed by Ginga in the 9.1-30 keV energy 
band (Iwasawa et al. 1993). r is actually only an upper 
limit to the variability timescales of the primary contin- 
uum. If the absorbing matter is spherically distributed, 
such a relatively rapid variability must be related only 
to the fraction of nuclear photons, which are transmit- 
ted without being scattered. For log(NH,tr) = 24.5, more 
than 2/3 of the 20-200 keV incoming photons are scat- 
tered (Matt et al. 1999b). This implies that the observed 
variability could be substantially diluted. If we suppose 
that the light curve in Fig. ^ is given by the sum of a 
constant plateau (scattered photons) and a variable con- 
tribution (unscattered photons), the latter has 2.6 times 
higher dynamical range, or, conversely, a r lower by this 
amount, r ~ lO'* s is still not exceptional among Seyfert 
galaxies. It is interesting to notice that NGC 4945 is one 
of the nearby galaxies with the smallest estimate of the 
central dark object mass (Mdo ^ 1-6 x IO^Mq; Greenhill 
et al. 1997). This suggests that NGC 4945 might be an 
absorbed version of NGC 4051 (Lawrence et al. 1985), ac- 
creting at a sizeable fraction (~ 10~^) of the Eddington 
rate. Alternatively, our estimate of the intrinsic nuclear 



variability dynamical range could be lowered if the cover- 
ing fraction of the absorber is much less than unity. The 
lack of any significant detection of Compton refiection in 
the broadband X-ray spectrum supports this hypothesis. 

It is remarkable that no spectral variability is associ- 
ated with these flux changes. This rules out that changes 
in the interposing absorbing medium are responsible for 
the observed flux variations. On the other hand, this 
suggests that the shape of the primary continuum in 
NGC 4945 is not significantly dependent on the inten- 
sity state. The infrared luminosity, (3 x 10'*'^ erg s~^ if 
calculated according to Mulchaey et al. 1994), and our 
extrapolation of the 2-10 keV nuclear intrinsic luminosity 
(3 X 10^^ erg s~^) lie well on the low luminosity end of 
the correlation observed in Seyferts for these two quanti- 
ties (Mulchaey et al. 1994). If Li^io kcv/Lboi 0.05, as 
typical for quasars with Li_io koV < lO'"^ erg s'^ (Elvis 
et al. 1994), the bolometric luminosity associated to the 
AGN is of the same order of magnitude as the observed 
infrared luminosity, and hence the AGN dominates the en- 
ergy output, at variance with the deductions of Genzel et 
al. (1998) on the basis of the ISO spectrum (see Marconi 
et al. 1999 for a possible explanation of this discrepancy). 

A high-energy cut-off in the primary nuclear contin- 
uum is measured for the first time. Its presence is sig- 
nificantly required in all the models adopted to describe 
the broadband emission, and in the PDS data alone as 
well. The cut-off energy cannot be, however, very well 
constrained, and values in the range 100-300 keV are 
possible. A reanalysis of the 50-500 keV OSSE spec- 
trum yields a steeper index than measured by BeppoSAX 
(Posse = 2.2±q;|), and also the BATSE average spectrum 
shows evidence for a high-energy cutoff above 100 keV 
(A.Malizia, private communication). Such cut-off ener- 
gies are well consistent with those measured in several 
Seyfert Is so far, both by OSSE (Zdziarski et al. 1995; 
Madjeski et al. 1995) and BeppoSAX (Piro et al. 1998; 
Guainazzi et al. 1999b; Guainazzi et al. 1999c; Perola et 
al. 1999). Also the intrinsic spectral index of the primary 
nuclear continuum (ranging between 1.4 and 1.7 in the var- 
ious models) is slightly but not exceptionally flat among 
Seyfert galaxies (Nandra et al. 1997; Turner et al. 1997). 

5.2. The 0.1-6 keV extended emission: unresolved 
discrete sources or truly diffuse emission? 

The total suppression of the transmitted component below 
~8 keV allows us to observe other spectral components, 
which are at least in part associated with the host galaxy. 
The emission in the 0.1-6 keV energy band is clearly ex- 
tended along the plane of the host galaxy, thus dismissing 
the possibility that the bulk of the intermediate X-ray 
emission is due to scattering of the nuclear radiation only 
(Iwasawa et al. 1993). In the MEGS image at least two 
sources are detected with a signal-to-noise ratio higher 
than 3. The spectral properties and luminosity of the 
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brightest arc consistent with a binary system accreting at 
sizeable fraction of the Eddington luminosity, if it is indeed 
associated with NGC 4945. Other explanations are, how- 
ever, not excluded by our data. If the contribution of these 
sources is taken into account, a significant extended emis- 
sion remains within at least the innermost 5 kpc, whose 
extent is almost symmetric with respect to the nucleus. 
This would lead to the conclusion that it is produced in 
a truly diffuse interstellar medium. There is, however, no 
way of determining the amount of the contribution of any 
unresolved underlying sources. The ROSAT/PSPC images 
of NGC 4945 suggest that at least 50% of the 0.1-2.4 keV 
photons arc produced by discrete sources. However, the 
extrapolation of the contribution of these sources in the 
MECS band yields a fraction of the observed flux varying 
between 15% and 50%, if a power-law with T comprised 
between 1 and 2 is assumed. The study of a large sample 
of early type galaxies with ASCA and BeppoSAX sug- 
gests that the bulk of the emission above 1 keV probably 
originates as the integrated emission from X-ray binaries 
(Matsumoto et al. 1997; Trincheri et al. 1999). A com- 
plete characterization of the properties of this extended 
emission should await the superior spatial resolution and 
sensitivity available with the scientific payloads onboard 
Chandra and XMM. What BeppoSAX data can clearly 
show is that both the iron line and the continuum above 
7 keV (where the presence of the transmitted primary nu- 
clear continuum starts to be dominant) are, by contrast, 
produced by a point-like region (at least as seen by the 
MECS). This strongly supports a nuclear origin for both 
components. 

5.3. The broadband X-ray spectrum 

The 0.1-10 keV continuum can be fit with the superposi- 
tion of a single temperature thermal emission from a colli- 
sionally excited plasma with kT ~ 3 keV and a warm scat- 
tered nuclear power-law, with a scattering fraction ~0.4%. 
This model has already been successfully employed to de- 
scribe the moderate resolution intermediate X-ray spec- 
trum of several Compton-thick Seyfert 2 galaxies (Turner 
et al. 1997; Guainazzi et al. 1999a; Vignati et al. 1999). 
Two main differences arc, however, present. The warm 
scattered power-law is seen through a substantial absorb- 
ing column density of ~3 x 10^^ cm~^. This might imply 
that the nuclear region is encompassed by absorbing mat- 
ter also on scales with are much larger than the dimen- 
sion of the Compton-thick structures, which almost sup- 
press the direct view of the nucleus. It is straightforward 
to associate it with the prominent dust lanes crossing the 
galaxy plane, although the presence of a starburst ring on 
a 100-pc scale (Marconi et al. 1999) provides a possible 
alternative source of absorption (Fabian et al. 1999). Re- 
cent HST NICMOS observations suggest that the active 
nucleus in NGC 4945 might be obscured along all lines 
of sight by matter with Nh ~ 10^-^-10^^ cm~^ (Marconi 



et al. 1999). On the other hand, the thermal component 
has a much higher temperature than typically associated 
with the X-ray emitting nuclear starburst in the Seyfert 2s 
(~500 eV). Again, this may simply be telling us that the 
single temperature description of this component is too a 
crude approximation. Scenario a), however, can hardly be 
reconciled with the fact that the bulk of the emission be- 
low 6 keV is extended on scales larger than a few kpc. The 
warm scattered power-law, in fact, dominates the fit in the 
whole 3-10 keV range, therefore also where the emission 
is extended. This is almost at odds with its supposed nu- 
clear origin. The power-law could of course represent only 
a phenomenological description of a different source of ra- 
diation, e.g. the contribution of a population of unresolved 
discrete sources. 

A formally equivalent description of the 
soft/intermediate X-ray spectrum is provided in the 
framework of the two thermal components template, 
used to fit the integrated bulge emission of X-ray faint 
early-type galaxies (Bregman et al. 1995; Kim et al. 1996; 
Matsumoto et al. 1997; Irwin & Bregman 1999). The 
temperature obtained with BeppoSAX for the hotter 
component is remarkably consistent with that measured 
on other early type galaxies, whereas the typical temper- 
ature of the soft component is slightly lower than usual 
(Matsumoto et al. 1997; Irwin & Bregman 1999), although 
consistent with that measured by ASCA in at least one 
case (NGC 4392; Matsumoto et al. 1997). The fact that at 
least 50% of the soft X-rays originate in discrete sources 
(see Sect. 3.2) is in agreement with the results on M 31 
(Primini et al. 1993; Irwin & Bregman 1999; Trincheri et 
al. 1999). On the other hand, Brandt et al. (1996) report 
a diffuse emission using the ROSAT/HRI on scales of 5' 
along the plane of NGC 4945, which is consistent with our 
reanalysis of the same data. Obviously, there is no way in 
this scenario to produce a neutral iron line with >1 keV 
EW. Actually, such lines have never been observed in the 
integrated spectra of early type galaxies (Matsushita et 
al. 1994; Matsumoto et al. 1997). The origin of the fine 
must therefore be connected with the nuclear activity. 

Recent studies suggest a possible alternative interpre- 
tation. In M 82 and NGC 253 (Cappi et al. 1999b), the in- 
termediate X-rays arc dominated by hot gas with kT ^6 - 
9 keV, probably associated with starburst superwind out- 
flows protruding from the galaxy disks. The observed iron 
abundances are strongly sub-solar. If this scenario is vi- 
able also for NGC 4945, these abundances can explain the 
lack of detection of fluorescent emission from highly ion- 
ized iron. The presence of a wind-blown cavity (Moorwood 
et al. 1996) strongly supports the existence of a starburst 
superwind in NGC 4945. Again, high resolution imaging 
with Chandra and XMM will provide invaluable contribu- 
tions to resolve this issue. 
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5.4- The iron line: signature of the nuclear absorber? 

The spectral fits do not require tlie presence of any com- 
ponent associated with Compton-reflection from the inner 
side of the putative molecular torus surrounding the nu- 
clear region, as observed in NGC 1068 (Matt et al. 1997) 
or Circinus Galaxy (Matt et al. 1999a). This rules out the 
most straightforward interpretation for the origin of the 
prominent (EW ^ 1 keV) fluorescent iron line. The 
centroid energy and intrinsic narrowness imply an origin 
from neutral or mildly ionized iron (Fe<xx). This is not 
in principle inconsistent with the line being produced in 
the same lukewarm medium responsible, in scenario a), 
for the almost energy-independent scattering of the pri- 
mary nuclear radiation. However, given also the problems 
of the scenario a), an appealing alternative explanation 
is that the line is produced in transmission by the same 
thick absorbing medium covering the nucleus (Leahy et 
al. 1989). Matt et al. (2000) compare the observed EW in 
the whole sample of the Compton-thick Seyfert 2 galaxies 
observed by BeppoSAX with the values expected in the 
pure transmission scenario, as calculated by the Matt et 
al. (1999b). In NGC4945, the observed EW is consistent 
with the model predictions, suggesting that the iron line 
is indeed due to transmission by the same absorbing mat- 
ter covering the nucleus. It is interesting to compare this 
outcome with that for the Circinus Galaxy, where the ob- 
served value (~60 eV against the total continuum; Matt 
et al. 2000) is more than one order of magnitude higher 
than expected from the column density of the absorbing 
matter (4 x 10^^ cm~^). In Circinus Galaxy a Compton 
reflection continuum is required to fit the broadband X- 
ray spectrum (Matt et al. 1999a), most likely providing 
the dominant contribution to the iron line. The difference 
between the NGC4945 and Circinus Galaxy case may be 
due to the geometry of the absorbing matter, allowing in 
the latter the direct view of the farthest side of the molecu- 
lar torus encompassing the nucleus. The reader is referred 
to Matt et al. (2000) for a more extensive discussion of 
this issue. 

6. Conclusions 

The main results of this paper can be summarized as fol- 
lows: 

a) the 0.1-5 kcV image profile is extended along the plane 
of the host galaxy. The extension is symmetric within 
the innermost 5', but at larger radii only the SE MECS 
profile shows a broader tail. This asymmetry is due to 
two discrete galactic sources which are present both in 
the soft X-ray ROSAT/PSPC and in the intermediate 
X-ray MECS/BeppoSAX fields. The profile perpen- 
dicular to the galaxy plane is possibly also extended 
within ±3', although this result is made more uncer- 
tain by possible residual inaccuracies in the BeppoSAX 
attitude reconstruction 



b) we measure the best high-energy {i.e.: >10 keV) spec- 
trum so far. Along with the lack of significant exten- 
sion in the MECS emission above 7 keV, the spec- 
tral fitting confirms the idea that the emission above 
10 keV is dominated by a nuclear non-thermal contin- 
uum, seen through an absorbing screen of Nh — a few 
10^"' cm^^. A flux variability with a dynamical range 
~60% and extrapolated doubling/halving timescales 
< a few 10^ s (similar to that discovered by Iwasawa 
et al. 1993 with Ginga) is not associated with any spec- 
tral changes, suggesting that it echoes a, possibly more 
extreme, variability of the intrinsic nuclear continuum 

c) from our estimate of the intrinsic nuclear power, we 
conclude that the energy output in NGC 4945 is likely 
to be dominated by the AGN, not by the starburst as 
suggested by infrared diagnostic (Genzel et al. 1998). 
NGC 4945 is therefore another indication that ab- 
sorbed AGN (and, possibly, also the so far elusive 
type 2 quasars) should be better searched in X-rays 
rather in the optical/infrared band 

d) the prominent iron line is consistent with fluores- 
cence from neutral or mildly ionized iron (<Fexx). 
There is no evidence for blending or intrinsic broaden- 
ing at a level higher than 200 eV. The bulk of the line 
(along with the continuum above 7 keV) is produced 
within about 1', strongly supporting a nuclear origin. 
The most likely explanation is fluorescence in the same 
absorbing medium, responsible for the nuclear absorp- 
tion, provided its covering fraction is lower than one 

e) the broadband 0.1-200 keV spectrum can be fitted 
with several statistically equivalent models. However, 
the most consistent explanation is in the framework 
of the two thermal components template, already em- 
ployed to fit the integrated bulge emission of early type 
galaxies, where the bulk of the soft X-rays and a size- 
able fraction of the intermediate X-rays are believed 
to be produced by unresolved discrete sources (mainly 
binaries). Alternatively, the extended emission above 
1 keV might be associated with a hot superwind-driven 
outfiow, as already observed in nearby starburst galax- 
ies. 
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